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1 introduction 



Gravitational waves in General Relativity are considered as perturbations of the space-time 
geometry which propagate with the speed of light. Contrarily to the behaviour of electro- 
magnetic waves in flat space, when gravitational waves pass through each other there is a 
non-linear interaction. Of course, such a difference arise from the fact that, unlike Maxwell's 
equations, Einstein's equations are highly non-linear. The situation is very different in the 
linearized theory of gravity where gravitational waves superpose without interaction. Even 
in the full Einstein-Maxwell theory. Maxwell's equations remain linear indicating non direct 
electromagnetic interaction between two pure electromagnetic waves. In that case, however, 
there is a non-linear interaction of the waves through the gravitational field generated by 
their electromagnetic energy. Although very small, the magnitude of such an interaction is 
similar to the magnitude of the interaction between pure gravitational waves. 

In order to find exact solutions of Einstein's equations describing plane wave collisions 
it is necessary to introduce several simplifying assumptions. Perhaps, the most restrictive 
simplification would be the assumption that the interacting waves have plane symmetry. 
Indeed, exact gravitational plane waves are very simple time dependent plane symmetric 
solutions of Einstein's equations That may not seem a very restrictive assumption if 
we intuitively think of plane fronted waves as approximations to spherical waves at large 
distances from their sources. However, the global features of the two type of waves are very 
different. In fact, exact gravitational plane waves exhibit two main surprising global features, 
namely: i) the absence of a global Cauchy surface, which is a consequence of the focusing 
effect that the waves exert on null rays 0, ii) the presence of a coordinate singularity on 
some hypersurface behind the wave. Such a singularity may be physically understood in 
terms of caustics that are formed in the regions where the null rays are focused . For pure 
gravitational plane waves, a whole three-dimensional congruence of parallel null geodesies 
may be focused into a line. For pure electromagnetic plane waves, the region of focusing 
may be a simple event. From this property it follows that we can conveniently consider the 
inverse of the focusing time as a measure of the strength of the wave. In particular, for 
a pure electromagnetic wave such inverse time equals the electromagnetic energy per unit 
surface of the wave. 

Beyond the assumption that the wave front is plane, the imposition of plane symmetry 
also requires that the magnitude of the wave is constant all over the entire plane. Further- 
more, if we only consider head on collisions, then it makes sense to impose the condition 
of plane symmetry globally. Actually, it may seem that there is not loss of generality if we 
concentrate only on head on collisions. The reasons are that the particular case of plane 
waves propagating in the same direction is trivial since the waves do not interact, they sim- 
ply superpose linearly |^ |^ . Also, we may think that it is always possible to make a Lorentz 
transformation in order to include oblique collisions. However, the assumption of global 
plane symmetry would severely blur the physical interpretation of the solutions obtained in 
this particular way. 

Nevertheless, besides their simplicity, one expects that exact plane waves may be relevant 
for the study of the strong time dependent gravitational fields that may be produced in the 
collision of black holes |^ or to represent travelling waves on strongly gravitating cosmic 
strings 0. In recent years these waves have been used in classical general relativity to test 
some conjectures on the stability of Cauchy horizons |]^, |T^, and in string theory to test 
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classical and quantum string behaviour in strong gravitational fields [TT|, |T2l |T3|. Their 



interest also stems from the fact that plane waves are a subclass of exact classical solutions 
to string theory |I5|, |1^ . 

It is known that when waves are coupled to quantum fields there is neither vacuum po- 
larization nor the spontaneous creation of particles. In that sense they behave very much 
as electromagnetic or Yang-Mills plane waves in flat space-time |]17| , |T^. Still the classical 
focusing of geodesies has a quantum counterpart: when quantum particles are present the 
quantum field stress-energy tensor between scattering states is unbounded at the Cauchy 
horizon, i.e. where classical test particles focus after colliding with the plane wave |^ . 
This suggests that the Cauchy horizon of plane waves may be unstable under the presence 
of quantum particles. The classical instability of the null Cauchy horizons of plane waves 
is manifest when non-linear plane symmetric gravitational radiation collides with the back- 
ground wave, i.e. when two plane waves collide. In this case the focusing effect of each wave 
distorts the causal structure of the space-time near the previous null horizons and either 
a spacelike curvature singularity or a new regular Killing-Cauchy horizon is formed. How- 
ever, it is generally believed that the Killing-Cauchy horizons of the colliding plane wave 
space-times are unstable in the sense that "generic" perturbations will transform them into 
spacelike curvature singularities. In fact, this has been proved under general plane sym- 
metric perturbations ||2^. Also exact colliding plane wave solutions with classical fields are 
known that have spacelike curvature singularities and which reduce, in the vacuum limit, to 



colliding plane wave solutions with a regular Killing-Cauchy horizon |21 



Note that the singularities derived from plane wave collisions are not the result of the 
collapse of matter but the result of the non-linear effects of pure gravity. Note also that the 
presence of a Killing-Cauchy horizon in a colliding plane wave space-time implies a break- 
down of predictability since the geometry beyond the horizon is not uniquely determined 
by the initial data posed by the incoming colliding waves. In fact, these extensions exhibit 
a wide variety, they may or may not contain space-like curvature singularities. Even they 
may contain time-like singularities, which can be avoided for certain observers in a space-like 
motion (see for a review on the subject). 

The interaction of quantum fields with colliding plane waves was first considered by 
Yurtsever ||2^ for the singular Khan-Penrose solution ^4], which describes the collision 
of two plane impulsive gravitational waves. In that case, an unambiguous "out" vacuum 
state was possible to define in a relatively simple way. More recently, Dorca and Verdaguer 
P3| , noticed that the presence of a Killing-Cauchy horizon in a non singular colliding 
plane wave space-time could be used to define an unambiguous "out" vacuum state related 
to the preferred Hadamard state introduced by Kay and Wald in more generic space-times 
with Killing-Cauchy horizons With this premise, Dorca and Verdaguer studied the 

interaction of quantum fields in a particular non-singular colliding plane wave space-time, 
the interaction region of which was isometric to a region inside the event horizon of a 

Later on the same premise was applied by Feinstein and 
which represents the head on collision of two 



Schwarzschild black hole [28, 29 



Sebastian 301 to the Bel-Szekeres solution 31 



electromagnetic plane waves with an interaction region isometric to the Bertotti-Robinson 
universe filled with an uniform electric field. In all these examples it was found that 
the initial state, defined to be the vacuum state in the fiat region before the arrival of the 
waves, contained a spectrum of "out" particles consistent, in the long wavelength limit, with 
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a thermal spectrum with a temperature inversely proportional to the focusing time of the 
waves. 

A further step in the study of the interaction of quantum fields with colliding plane 
waves is the computation of the expectation value of the stress-energy tensor. However, the 
non perturbative evaluation of the stress-energy tensor of a quantum field in a dynamically 
evolving space-time is generally a difficult task. Even when the exact modes of the quantum 
field equation are known it may not be possible to perform the mode sums in order to get the 
quantum field two point function or, more precisely, the Hadamard function, which is the key 
ingredient in the evaluation of the stress-energy tensor. This problem was first considered by 



Yurtsever |23] for the Khan-Penrose solution [24|. In that case it was possible to determine 



the behaviour of the stress-energy tensor near the singularity of the interaction region. It 
was shown that for the conformal coupling case (i.e. = 1/6) the energy density and two of 
the principal pressures were positive and unbounded towards the singularity. This problem 
has been also considered by Dorca and Verdaguer in the mentioned above non-singular 
colliding plane wave spacetime with an interaction region isometric to an interior region 
of a Schwarzschild black hole. As in the case of Yurtsever for the Khan-Penrose solution, 
the expectation value of the stress-energy tensor was calculated in the state representing 
the Minkowski vacuum in the flat region before the arrival of the waves. This value was 
first computed in a region close to both the Killing- Cauchy horizon and the topological 
singularities, the folding singularities, that the colliding plane wave space-time contains ^ 



In that particular region, the calculations were simphfied due to the blueshift effect on the 
energy of the initial quantum modes as they reached the Killing- Cauchy horizon |33]. In 



a recent work [0], an approximation procedure was proposed by the author in order to 



calculate such an expectation value throughout the causal past region of the collision center. 
That approximation has been also applied by the author in the Bel-Szekeres space-time |p5|| . 
In all of these calculations, it was found that the stress-energy diverged as the Killing- Cauchy 
horizon was approached. The rest energy density was positive and unbounded towards the 
horizon. Two of the principal pressures were negative and of the same order of magnitude 
of the energy density. It was also pointed out that such a behaviour suggested that the 
non singular Killing- Cauchy horizon is indeed unstable under quantum perturbations and a 
curvature singularity would be the general outcome of a generic plane wave space-time when 
backreaction is taking into account. Note that this is a non perturbative effect, it is the 
result of the nonlinearity of gravity, since gravitational waves in the linear approximation 
do not polarize the vacuum. In fact the vacuum stress-energy tensor of a quantum field in 
a weakly inhomogeneous background was computed by Horowitz [Q, and it is easy to see 
that such tensor can be written in terms of the linearized Einstein tensor only [^, which 
vanishes for gravitational waves [p8| , 

The plan of the paper is the following. In section 2 the geometry of plane fronted waves 
and colliding plane waves is briefly reviewed. In section 3 the mode solutions of the scalar 
field equation are given throughout the causal past of the collision center in the four different 
regions of the space-time. It is explicitly seen that an exact expression for these mode 
solutions can be easily found everywhere but in the interaction region. However, following 
the directions of two previous works [^, |35|, we will give an adequate approximation for 
the field solutions in the interaction region. Such an approximation will allow us to recover 
the singularity pattern of a Hadamard function, which is essential in order to perform a 
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correct regularization. In section 4 a brief summary of the main point- splitting regularization 
formulae is given. Finally, in section 5, a procedure to calculate the Hadamard function is 
introduced. Also, are given the basic details to compute the vacuum expectation value 
of the stress-energy tensor in the causal region of the collision center. A summary and 
some consequences of the results are discussed in section 6. Also, some bitensor covariant 
expansions have been stored in the Appendices. 

2 Description of the geometry 

In this section we start with a review of the geometrical properties of gravitational plane 
waves, and we follow with a geometrical analysis of the head on collision of two gravitational 
plane waves. 

2.1 Gravitational and electromagnetic plane waves 

The appropriate notation to describe plane waves is that of Newman and Penrose (1962) 
|[40|| . We introduce a thetrad on null vectors, namely, two real vectors ra^, and a complex 
null vector and its complex conjugate rh'^, such that ra'^Z^ = 1, m^ffi^ = —1, and they 
satisfy the completeness relation, 

gtiu = liin,y + n^ly - m^rhy - rh^my. 

Then we can write the ten independent components of the Ricci and the ten independent 
components of the Weyl tensor in this thetrad basis. This turns to be physically convenient 
since the components of the Weyl tensor in the thetrad basis have a well established meaning 
of Coulomb components or transverse and longitudinal wave components in the null direc- 
tions n^, |¥T]. This interpretation is also very useful to analyze the interaction of two 



of such waves since it is possible to align the basis vectors n^^ and with the propagation 
directions of the two waves. Furthermore, if we consider the interaction between transverse 
waves, the problem reduces to find the interaction between the respective transverse Weyl 
components. 

In the context of plane waves we introduce the class of pp-waves, i.e plane fronted grav- 
itational waves with parallel rays. They are defined by the property that they admit a 
covariantly constant null vector field, which is possible to interpret as the rays of the wave 
and it may be identified with the thetrad vector l^. There exist a family of 2-surfaces or- 



thogonal to that may be interpreted as wave surfaces [0. It is convenient to use a null 
coordinate u, defined by = u^^, and the metric can be written in the Kerr- Shild form as, 

ds^ = 2 dudr + H{u, X, Y) du^ - dX^ - dY^, (1) 

where the coordinates X and Y span the wave surfaces. These type of metrics are either of 
algebraic type N, or conformally flat and the only non-zero component of the Weyl tensor, 
in the thetrad basis, is the transverse wave component in the direction. The modulus 
term of this component can be conveniently interpreted as the amplitude of the wave and 
its phase term as the polarization of the wave. 
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If we suppose that these space-times are solutions of Einstein's equations, the function 
H in the metric (|I]) must satisfy the differential equation. 



H^xx + Hyy = 0. 

Recall that this is a linear differential equation which means that distinct solutions of pp- 
waves (|I]) may be superposed, i.e. two independent pp-waves, propagating in the same 
direction do not interact 0, ^ . Another remarkable property of pp-waves, which was pointed 
out by Yurtsever |Q, states that any gravitational wave space-time that is flat before the 



arrival of the wave and returns to a perfect flatness after the wave passes, i.e. a sandwich 
wave, it is necessarily a pp-wave. 

In Einstein-Maxwell theory, the particular class of plane-symmetric waves are defined to 
be pp-waves in which the curvature field components are the same at every point of the 
wave surfaces, i.e. the transversal X, F-planes. There is only one non-null component of 
the Ricci tensor and, as mentioned above, only one non-null component of the Weyl tensor. 
Both components depend only on second derivatives of the function H{u, X, Y) with respect 
to the transversal coordinates X and Y. Therefore, the plane symmetry condition implies 
that H{u, X,Y) is quadratic in X and Y, i.e.. 



H{u, X, Y) = hniu)X^ + 2 hi2iu)XY + h22iu)Y^. 

In that case the non null components of the Ricci tensor, in the thetrad basis, is given by 
the combination {hu + /i22)/2 and the non null component of the Weyl tensor, in the thetrad 
basis, i.e. the transversal wave component in the direction is given by the combination, 
I (/ill — + i2 hi2). Then, depending on the values of hij{u), we will distinguish between 
two type of waves, namely, i) pure gravitational waves, when hu = —h22, and therefore the 
Ricci tensor is exactly zero and we have a vacuum solution of Einstein's equations. If in 
addition hu is proportional to hi2, then the wave has constant linear polarization, ii) pure 
electromagnetic waves when hu = h22 and hi2 = 0, and therefore the Weyl tensor is exactly 
zero. 

In order to consider the collision and interaction of plane waves it is first convenient to 
change to Rosen-like coordinates {u, v, x, y), where u and v are two null coordinates, x and 
y two transverse coordinates, and the line element (|l|) can be transformed into, 

ds'^ = 2 dudv - e"^ (e^ cosh Wdx"^ - 2 sinh Wdxdy + e"^ cosh Wdy^'^ , (2) 

where U, V and W are functions of u only. In the case of linear polarization, it is always 
possible to set = 0. 

However, the line element in Rosen form, always contains a coordinate singularity 
on some hypersurface behind the wave. This singularity appear as a consequence of the 
focusing effect which plane waves exert on null geodesies. Physically it may be interpreted 
as a caustic produced in the region of null ray focusing 0. Related to this focusing property 
of plane waves, as was first pointed out by Penrose 0, there is the fact that plane wave 
space-times do not contain global Cauchy surfaces, i.e it is not possible to set up initial 
values for a plane wave on any global space-like surface which lies entirely at the causal past 
of the wave front. 
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2.2 Colliding waves 

The colliding plane wave space-time consists of two approaching waves, regions II and III, 
in a flat background, region IV, and an interaction region, region I. The two waves move 
in the direction of two null coordinates u and v, and since they have plane symmetry in 
the direction of the transversal coordinates x and y, the interaction region retains a two- 
parameter symmetry group of motions generated by the Killing vectors and dy. At each 
point of the interaction region, the null directions, along coordinates u and v, are orthogonal 
to the planes spanned by and dy. We may align the thetrad vectors n'^ and l'^ with these 
two null directions in such a way that, = u^^/L and = v^^/N, being L and two 
functions which do not depend on coordinates x and y. 

The general line element for the interaction region can be written in Szekeres form as 

ds'^ = 2e-^dudv - e"^ (e^ cosh W + 2 sinh Wdxdy + e"^ cosh Wdy^) , (3) 

where M, U, V and W are functions of u and v only. These functions are constrained by 
Einstein's equations together with a suitable set of initial conditions at the boundaries with 
the plane wave regions II and III. It is possible to set in general, 

e-^ = a{u) + b{v), (4) 

with a{u) and b{v), two arbitrary monotonically decreasing ivMciions that may be considered 
piecewise C^. Recall that the line element (0) is rather convenient because it is similar to 
that for a single plane wave in Rosen form @, except that the functions M, f/, V and W 
are functions of both coordinates u and v in region I, of coordinate u alone in region II and 
of coordinate v alone in region III. In region IV all these functions are constant and may be 
set to be zero by an appropriate coordinate parametrization. 

The appropriate junction conditions across the null hypersurfaces m = and f = 0, are 
those of O'Brian and Synge P^. These conditions demand that only the components, g^^, 



g'^^gijfl, g^^gijfl, should be continuous through the null hypersurfaces. Notice that O'Brian 
and Synge conditions are weaker than the usual Lichnerowicz conditions, which require that 
the metric tensor is and piecewise at least C^. In particular, Lichnerowicz conditions 
exclude the possibility of impulsive gravitational waves. 

The function U is appropriately set by the boundary conditions (see for instance P^): 

a{u) = — 0; '^(0) — 0; 

h{v) = ^; < 0, and 6(0) = 0. 

Thus, the values of exp(— f/) are 1 in region IV, l/2 + h{v) in region III, l/2 + a{u) in region 
II and a{u) + h{v) in region I. The statement that U is an smooth function across the null 
boundaries requires that the function a{u) must have the same form in both regions I and 
II. Similarly, the function h{v) must have the same form in regions I and III. Furthermore, 
being a[u) and h{v) monotonically decreasing functions, it is always possible to express them 
in the form, 

a{u) = \-Q{u) (ci ur , h{v) = \-Q{v) (C2 vr , (5) 
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for a certain positive values of the parameters Ci, C2, ni, ^2 and where 0(C) is the usual 
step function. Recall that this expressions for a{u) and b{v) apply globally. Observe that, 
since functions a{u) and b{v) are monotonically decreasing, it is inevitable that a singularity, 
either a curvature singularity or a Killing Cauchy horizon, at a{u) +b{v) = is formed. The 
two constants Ci and C2 in are a measure of the strength of the waves. In fact, the larger 
these parameters are, the sooner the singularity is produced. It is also possible to use further 
rescaling on null coordinates u and v to set ci = C2 = 1, and with this rescaling, as depicted 
in Fig. 1, we have a symmetric wave collision. Once a rescaling is used, to set ci = C2 = 1, 
then a measure of the strength of the waves is absorbed into the function exp(— M) in line 
element @. In fact, exp(— M) will be proportional to L1L2, where the two new parameters 
Li and L2 are directly related to Ci and C2 as Li = cf^ and L2 = C2 ^• 

Once a value for function U is set, Einstein's equations together with suitable initial 
conditions at the boundaries with the null hypersurfaces u = and v = 0, will determine 
the values for the remaining functions, i.e V, M and W. For the particular case of linearly 
polarized plane waves, we may set W = 0. In that case, one possible family of solution of 
Einstein's equations which satisfy the boundary conditions is the Szekeres Class of solutions 



451, for which. 



^viu,v) ^^^^ ^y,+p, I 1^12 + 5 + ^1/2 _ a ) fyi/2 + a + ^1/2 - h\ , (6) 



. ffl + ftV(Pi+P2)2-l]/2 / \ -2piP2 



(l/2 + a)?'i/2(l/2 + 6)Pi/ 



(7) 

where pi and p2 are two real parameters related to the parameters ni and n2 in the definition 
of functions a{u) and h{v) (H) by. 



p =2 1-- , p =2 1- 

\ riiJ \ n2 

In order to satisfy the boundary conditions, it is necessary that, > 2 with i = 1,2, and 
therefore the constants pi are restricted to the range of values 1 < pj < 2. For instance, the 
particular values pi = P2 = ^ correspond to the Khan- Penrose solution In this family of 
solutions, there is always a scalar polynomial curvature singularity in the interaction region 
on the surface a{u) + b{v) = 0. 

We will restrict in the following sections to the Szekeres Class of solutions. The colliding 
plane wave space-times can be briefly described as follows (see Fig. 1): The general space- 
time contains four regions, namely: two single plane wave regions (regions II and III) moving 
in a Minkowski fiat background (region IV) and an interaction region (region I). These four 
space-time regions are separated by the two null wave fronts u = and v = 0. Namely, 
the boundary between regions I and II is {0 < m < 1, v = 0}, the boundary between 
regions I and III is {u = 0, < v < 1}, and the boundary between regions II and III with 
region IV is {u < 0, v = 0} U {u = 0, v < 0}. The interaction region is bounded by 
a scalar polynomial curvature singularity at the surface a{u) + b{v) = 0. Region I meets 
region IV only at the surface u = v = 0, and the plane wave region II or III meets the 
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singularity of region 1 at V = {u = 1, v = 0} or V = {u = 0, v = 1}, respectively. 
Observe that, the single plane wave regions 11 and 111 always contain a singularity at the 
hypersurface u = 1 for region 11 and v = 1 for region 111. For the particular case of the Khan 
and Penrose solution, which corresponds to the values rii = ^2 = 2, these singularities are 
just coordinate singularities. For all the other cases, i.e. ni,n2 > 2, these singularities are 



non-scalar curvature singularities ^g]. However, in what follows we will refer generically 
to them as folding singularities. The reason for this terminology comes from the behaviour 
of null geodesies near them. Namely, only a null measure set of null geodesies which enter 
into the plane wave regions 11 or 111 hit the singularities, all the rest are folded to end up 
into the interaction region. Such a behaviour suggest that the whole singularity at m = 1 in 
region 11 (or v = 1 in region 111) should be identified with the surface V = {u = 1, v = 0} 
(or V = {u = 0, V = 1}) in region 1 (see for instance [25] for a 3-dimensional plot of a 
space-time of this type). 



3 Mode propagation 

For simplicity we will consider in this section a massless scalar field, which satisfies the usual 
Klein-Gordon equation, 

□0 = 0. (9) 

Following the directions of the approximation procedure introduced in the previous works 
34| , |35| , we will be interested in the value of the quantum field (j) all over the causal past 



region of the collision center. The reason is essentially because the calculations can be 
greatly simplified in this region. We will start with the field solution in the flat region prior 
to the arrival of the waves, which is chosen to be the usual vacuum in Minkowski space- 
time. This vacuum solution will set a well posed initial value problem on the null boundary 
E = {m = 0, V < 0} U {u < 0, w = 0}, by means of which a unique solution for the field 
equation can be found throughout the space-time, i.e., in the plane wave regions (regions 
II and 111), and in the interaction region (region 1). However, although it is rather easy to 
find the solution of the field equation in regions 11 and 111 which matches smoothly with the 
boundary conditions, it turns out to be a difficult problem for the interaction region. The 
reason is essentially due to the intrinsic differences between the geometry of the plane wave 
regions and the interaction region. In fact, as mentioned in section 2, the plane wave regions 
are either conformally flat or type N in the Petrov classification, but the interaction region 
can be more generic. We will refer, from now on, to this problem as the mode propagation 
problem. 

We will consider the line element, 

ds^ = 2e-^^("'^)rfMrft; - e-^("'^) (e^^^'^^dx^ + e'^^^'^^ciy^^ , (10) 

which applies globally to the four space-time regions, and where the functions U, V and 
M, can be directly read off (|D-(0). Then, the field equation can be separated in a plane- 
wave form solution for the transversal coordinates x and y, with kx and ky, respectively, 
as separation constants. This plane- wave separation is just a trivial consequence of the 
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translational symmetry of the space-time on the planes spanned by the Kilhng vectors 

(11) 



and dy. The field solution is thus, 



V, X, y) = e^("'")/2 f^^^ ^) ^ik^x+ikyv^ 
where the function f{u,v) satisfies the following second order differential equation, 



,uv_ I ^ M+U 

e-^/2 ^ 2 



2„V 



(12) 



From now on, we will refer to Q{u,v) as the potential term. In the fiat region (region IV) 
this potential term is simply, 

njyiu,v) = L^L2ikl + kl), (13) 
where we have used that the functions U and V in (|10]) are zero in the fiat region, and the 



(14) 



function exp(— M) = 2L1L2. Using (pTSf), equation ( |12D can be solved as 



f{u,v) 



-i2k+u—i2k—v 



where k± are two new separation constants with dimensions of energy, and we define for 
convenience two dimensionless constants as k± = \fL^L^k±. These new separation constants 
are directly related to the previous ones k^ and ky by, 

Ak+k. = kl + kl. (15) 

The field solution in region IV reduces, thus, to the usual Minkowski plane wave solution, 
i.e, 



(t)k{u,v,x,y) 



2k_{2'Kf 



-i2k^u—i2k-v+ikxX+ikyy 



(16) 



These modes are well normalized on the null hypersurface E, which is the boundary of the 
plane wave regions II and III with the fiat region IV. Even though S is a null hypersurface, 
a well defined scalar product is given by (see [^] for details). 



—i / dxdy 



(j)ldu4>*2 



du + 



D = 



(f)idy(j)*2 



dv 



u=0 



(17) 



The modes ( p^Gf ) will determine on S a well posed set of boundary conditions for modes in 
regions II and III. There, the potential term in equation (|I2[) is simply. 



Q,{u,v) = ^e-''^+''^{kle-''^ + ky' 



where the label i = II or i = III in the functions U, V and M, stands for their particular 
values in the plane wave regions II or III. Then, the solution of equation (|T^ in regions II 
and III with the boundary conditions imposed by the fiat modes (|1^) on the hypersurface 
S, can be easily found as. 
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-i2k-v-iAn{u)/{2k- ) . 



in region II, 



(19) 



-i2k+u-iAiii{v)/{2k+) . ■ 



in region III, 



where the generic function Ai{(), with i =11, III is given by. 
Therefore, the well normalized "in" modes in regions II and III are. 



(20) 



COSM 
1 

( 

COSf 



~i2k^v-iA^^{u) /{2k- ) . 



-i2fe+u— iyljjj(f)/(2fc+). 



in region II, 



in region III, 



(21) 



Now the Cauchy problem is well posed on the boundary = {m = 0, < w < 1} U {0 < 
u < 1, V = 0} between plane wave regions II and III and interaction region. Notice that the 
initial modes (^1]) are well normalized on the boundary S between the flat region and the 
plane wave regions, and this means, from general grounds, that they remain well normalized 
on the boundary between the plane waves and the interaction region. This can be seen 
explicitly using a well defined scalar product on the hypersurface Si, which similarly to (|T7p 
is, 



—i dxdy 



-Uiii 



du + 



dv 



u=0 



(22) 

The correct normalization of modes ([2T| ) follows easily from the fact that functions Ai{x) in 
( PPD are unbounded at the folding singularities, V = {u = 1, v = 0} and V = {u = 0, v = 

!}• 

We have now to solve equation (|T^) in region I with the boundary conditions imposed by 
(pi]) on the lines Si= {^ = 0, 0<t'<l}U{0<M<l,t; = 0}, which are characteristic lines 
for the partial differential equation (|T^. Thus, the only independent boundary conditions, 
i.e. the Cauchy data, are the initial values of the function f{u,v) on them (the normal 
derivatives of the field on the characteristics, which are usually part of the Cauchy data, are 
determined by the values of the function f{u,v) itself, see for details). Furthermore, if 
we recall that the collision center is determined by the condition u = v, then the subset of 
Cauchy data that affects the neighborhood of the collision center, i.e its causal past, lies on 
the segments, Si = {u = 0, < v < s} Li {0 < u < s, v = 0} 0, where + s"^ = 1 and 
rti, n2 are defined in (^. We will denote the causal future region of these Cauchy data (or 
equivalently, the causal past region of the collision center) by region S, see Fig. 2. 

We can determine the behaviour of the solutions of equation ( P^ in the region near the 
singularity using the following change of coordinates. 



a 



b{v) + a{u), p = h{v) — a{u). 



(23) 
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Then equation (|T^) becomes, 



l^^-Up + nia,p)f = 0, (24) 

where the potential term is, 

ab ^ ^ 



n{a, p) 



4a2 



(25) 



and where we have used for convenience two dimensionless separation constants ki and A;2; 
directly related to k^ and ky by. 



fci = yLiL2 fcx, k2 = \lLiL2ky. (26) 

Observe that this equation has a singular point at a = due to the term cr"^ on the potential 
n((T, p). In fact, we have written the functions exp(— M) and exp(±V^) in the following way, 

e-*^ = a"e-*^ e^ = a^e^, a = (/J^ - l)/2, (3 = p^+p2, (27) 

as can be directly seen by (§|)-(|3)5 where the new functions exp(M) and exp(±V^), are 
bounded at cr = 0. It is straightforward to see from ( pT]) that a ± /? + 1 > 0. Thus, the 
singularity term in equation (|2^ ) is only given by the factor 0""^. 

Observe also that near the singularity a = a{u) + b{v) = 0, the potential term ( pSj ) blows 
up due to the factor but it is perfectly smooth on coordinate p = b{y) — a{u). Thus, 
the relevant features are related to coordinate a only. Therefore, we may not expect any 
physically remarkable difference if we substitute coordinate p in (|25|) by its constant value 
p = Ps = b{s) — a{s) when the singularity is approached from inside region S, i.e. when 
u, V ^ s such that s"'^ + s"^ = x. Then, we will expect a solution of equation (|2^) of the 
type /(cr, p) = g{a) explik^p), being k^ a separation constant and where the new function g 
depends only on coordinate a and satisfies the differential equation, 

g,aa + p(o-,Ps) + kl\ g = 0. 

We could now attempt to solve this Schrodinger-type of differential equation. For instance, 
we could try to find a solution by a WKB expansion. Unfortunately, the initial conditions 
for equation (p^ ) are given on the boundaries Si and there the coordinates a, p are not well 
defined. This means that, even though coordinates a and p would allow us to find a solution 
near the singularity a = 0, they cannot be used to find a solution near the boundaries with 
the plane waves, and therefore they are not convenient for the mode propagation problem. 
Nevertheless, it will be useful in what follows to know that this type of solution can be found 
near the singularity. 

In fact, we should return to equation (|T^) in coordinates u, v and we should recall that we 
are only interested in finding the field solution in the neighborhood of the collision's center. 
Thus, the potential term of equation (0) is restricted into region S, and the boundary 
conditions are given on the segments Sj (see Fig. 2). Then if, in analogy with the change 
(p3D , we use the notation a = , such potential term is simply given by, 

Q{u,v) = ^^ + ie-^V-^- {e-Hl + a'^eH^ , (28) 
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with a, (3 and e~^, defined in (pT]). Observe that functions e~^, e~^ and a{u)h{v) 
are perfectly smooth in the entire region S. However, ( pSD blows up when the term a = 
1 — u"'^ — v^^ goes to zero at the point u = v = s, where s"'^ + s"^ = 1. Taking all of this 
into account, it would be useful to find an adequate approximation for the potential ( p8l) 
such that: i) it preserves the essential features of (^), ii) it allows us to find a solution of 
equation (|12D which smoothly matches with the boundary conditions at Sj. 

In order to find such an approximation it would be rather convenient to introduce a 
new set of coordinates naturally adapted to this problem. Recall that coordinates (p3D are 
convenient to solve equation (|T2|) close to the singularity, but they are not well defined 
globally and particularly on the boundaries with the single plane waves. Nevertheless, since 
a = a{u) + h{y) behaves near the singularity in region 5 as, cr ~ {u — s) + bg {v — s), 
with ds = d(s), bg = b{s), we may define in analogy with ( p3D a new time-like and space-like 
coordinates as. 



t 



S) + CLs [U 



[a.bg 



-1/2 



a. [u 



{dsbs 



-1/2 



(29) 



where in fact, t and z apply globally and close to the singularity in region S they behave 
like coordinates (pSj). Indeed, t — *• a{dsbs)~^^'^ and z — > (p — ps){dsbs)~^^'^ as the singularity 
u = V = s is approached. The factor (dsbs)~^^^ is used to ensure that such a coordinate 
change represents a Lorentz transformation. 

We can easily find an expression for the function a{u, v) in coordinates t, z in the neigh- 
borhood of the singularity u = v = s hy expanding a{u, v) in terms of a Taylor series 
centered in m = w = s, with u — s and v — s expressed, from (p9[) , in terms of coordinates t 
and z. The result is, 



o{t,z) 



t (dsbs 



a/2 



+ E 



t - z 



+ 



rt + z 

if} 



{asbsT'^ 



(30) 



Observe, however, that in region S we expect the physical remarkable effects to occur near 
the singularity u = v = where the potential term ( PBD grows unbounded as a{t, z) 0. 
Recall that in the neighborhood of (j(t, 2;) = the relevant term in (|30| ) is the first order term 
t {dsbsY^"^, which is insensitive to coordinate z. Therefore, we may not expect any physically 
remarkable difference if instead of working with the function a{t, z), which depends on both 
coordinates t and z, we change to a new function a{t) depending only on coordinate t and 
directly derived from (|30| ) by taking z = throughout the entire region S, i.e.. 



ait)=tidsbsY^' + Yl 



asbsT'^ 



n=2 




(31) 



Observe also that the rest of the functions which appear in (pH]), i.e. e^^''"''"\ q-m{u,v) ^^^^ 
d{u)b{v), are perfectly smooth all over region S. Thus, we may not also expect any physically 
remarkable difference if we substitute these functions, which also depend on both coordinates 
t and z, by approximate functions depending on coordinate t alone. The procedure would be 
equivalent to the one used for a{u, v). Namely, given a smooth function ipi^u, v) in region S, 
we expand it by a Taylor series centered in the singularity u = v = s. Then, we substitute 
u — s, V — s hj coordinates t, z using the definition (p9|) and finally we take z = 0. The 
result would simply be. 
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' (32) 



2/ n! 



Recall that the boundary conditions for equation ([12|) restricted in region iS, lie on the 
segments Si and since the endpoints of these segments are not close to the singularities V 
or V\ we may suppose that the boundary conditions are not significantly different from the 
flat boundary conditions. This fact, together with the prior replacements of m, f -dependent 
functions by related t-dependent functions in the potential term (^Sf ), suggest that we may 
reproduce the main physical features if we replace the four-region mode propagation problem 
by a much simpler problem. Such a problem would consist in solving equation (|T2|) with 
a potential term (|28| ) depending on coordinate t alone and where the boundary conditions 
would be imposed by the flat Minkowski modes (|16D below the hypersurface {T = 0, —s< 
Z < s}, where we denote T = u + v,Z = v — u (see Fig. 3). Recall, however, that all this 
discussion is absolutely non applicable when a solution for equation ([T^) in the neighborhood 
of the folding singularities V or V' is required. This is non applicable because the boundary 
conditions (pTF) are unbounded as the folding singularities are approached and therefore these 



boundary conditions strongly differ from their counterparts (16) in the flat region. 

In order to solve this Schrodinger-type problem, rather than relying on the discussed 
approximations for the exact fleld equation ([l2|) , we will rewrite a new fleld equation using 
an adequate approximation for the line element throughout the entire causal past of the 
collision center. This new approach, which may seem redundant in the case of the particle 
production problem, is absolutely necessary when the renormalization of the stress-energy 
tensor is discussed. This is essentially because the process of renormalization involves the 
subtraction of the inflnite divergences that arise from the formal deflnition of the stress- 
energy tensor, and these divergences can be expressed as entirely geometric terms, which 
are independent of any possible approximations in the fleld equation. This means that in 
order to recover the geometric divergences in the stress-energy tensor, any approximation in 
the fleld equation must be related to a suitable approximation in the space-time geometry. 

The approximations just discussed above can be essentially recovered by changing the 
exact colliding plane wave line element throughout the causal region of the collision's center 



([lOD into a related line element. This new line element would basically consist in replacing 
the metric coefficients a{u,v) = e'^^"''^), a(M, t;)±^e±^("'^), a(M, t;)-"e-^^("'^) in (0), which 
depend on both coordinates u and v, by new functions of coordinate t alone in the sense 
described in (pl|)- (|3^) , i.e., 

dsf = i cr{tr e-^^W {dt' - dz') - a{tY+^ (e^^ dx' + ait)"^^ e'^^ dy^) . (33) 

We will suppose that the line element (|33D applies all over the causal past of the collision 
center, not only in the interaction region but also through the plane wave regions II and III 
in the sense of Fig. 3. The plane wave collision starts aX t = t^ = —s {bs + dg) (as6s)"^/^ but 
to avoid smoothness problems derived from such an approximation, we will suppose that 
( |33D applies exactly on a range t^ + e < t < 0, for a certain e > 0. In the range to < t < tQ + e, 
as described below, we will interpolate a line element which smoothly matches with the flat 
space at t = to- Nevertheless, the particular details of this matching will not affect the main 
physical features. 
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The exact field equation for this approximate space-time is, 



(□ + ei?)0 = O, (34) 

where it is necessary to consider a couphng curvature term in the field equation because, 
although the exact space-time is a vacuum solution, we have a bounded nonzero value for R 
in the approximated space-time. In order to solve this new field equation, we start rewriting 
the line element (^) in the following general way, 

ds' = ihhh) de' - dz- _ j^:^^ dx- _ ^Jjl]^ dy\ (35) 

where the /, are functions of coordinate t alone, which for values oi + t < t < can be 
straightforwardly determined by direct comparison with (^3]) as 



flit) = le-^W-*W a(t)(^-^)V2, /|(t) = ^_,m-m) o-(t)('3+l)^/^ fiit) = a\t), (36) 



where definitions { ^7\ ) have been used. For values t < to we take /i(t) = f2{t) = 
fsit) = 1, which correspond to their values in fiat space. Finally, in the interval to < t < 
to + e, we smoothly interpolate each fi(t) [i = 1,2,3) between these values. Also, in order 
to prevent singularities in the field equation, we conveniently reparametrize coordinate t, by 
t*{t), as follows. 

Now, we use the following ansatz for the field solutions, 

0;. = /^(r)e^^-^+^'=^^+^'=^^ (38) 

where the plane wave factor in coordinates x, y is related to the translational symmetry of 
the space-time along the transversal directions x, y, and the plane wave factor in coordinate 
z is just a consequence of our approximation. Then equation (p^ ) directly leads to the 
following Schrodinger-like differential equation for the function h{t*), 

V,. + u'it) h = 0, V{t) ^ co\t) = m + flit) kl + fiit) kl + flit) kl (39) 
where the function /o(i) stands for, 

flit) = m)f2it)fsit)] ^R. (40) 

Observe from (||) and (|27|), that the range of the parameter j3 in ( pG] ) is < 8. Then, we 
always have that either (/5 + 1)^/4 < 1 or (/? - 1)^/4 < 1. 

Now, it would be rather convenient to discuss whether or not it makes sense to solve the 
differential equation (|39D in terms of a WKB expansion. Observe that, when the parameter 
/3 = ±1 in (PBf), then the frequency u;(t) goes to a non-zero constant value towards the 



singularity a = 0. Then the short-wavelength condition, i.e. uj~^d/dt* Intu -C 1, holds and is 
particularly accurate near the singularity. Hence, equation (^) admits everywhere a WKB 
solution. 
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However, when (3 ^ ±1, the exponents of the cr-terms in the functions /j in (|36|) are 
all non-zero. Thus, all the functions /j vanish at the singularity a = 0. Also the function 
/q = ^ i? (/1/2/3) vanishes at a = 0. This can be easily seen because /o can be expressed in 
terms of /i, /2 and /a, as 



ll A- ll A- — - ( ll A- ll 4. ll\ 

fi f2 h 2 l/i h h 



and substituting the values (^) we have that /q is proportional to a. Therefore, the fre- 
quency term uj{t) in (|39D goes to zero towards the singularity a = 0. Hence, the short 
wavelength condition does not hold in the region close to the singularity and thus a WKB 
solution for equation ( p9| ) is not directly appropriate in this region. Nevertheless, we can 
still use a WKB solution, which at least would be appropriate if we are not close enough to 
the singularity cr = (recall that a = means from (pT]) that t* — » 00). Observe, however, 
that such a WKB expansion may still give information near the singularity. The reason lies 
upon the fact that we can interpret this WKB solution as a limit case of a set of WKB 
expansions which are truly appropriate near the singularity. Such a limit process may be 
briefly introduced as follows: we add a small quantity to V{t) in (^). Then, the short 
wavelength holds and it reduces to (dt/dt*) dV/dt <^ 2lj^, where now the frequency uj(t) is 
non-zero and dt/dt* = fslt) towards the singularity. Hence, the differential equation 
(|39D admits a WKB solution everywhere. We could then perform the involved calculations 
and finally set the small quantity that we have introduced in V{t) to zero. Mathematically, 
adding this small quantity to V(t) would be equivalent to subtract a small amount 6 > 
into the exponents of the a-terms of functions /i and /2 in (|36[) , i.e.. 



which is also equivalent to slightly change the line element (^) by including a factor a(t)^^^ 
to the metric coefficients gt*t*, 9zz only. Then, the function /q is finite at the singularity and 
is given by /g = 20^6^5 + 0{a). 

Finally, using a WKB expansion, the mode solutions (pk which reduce to the fiat mode 
solutions in the region prior to the arrival of the waves, are 

" 1/2 (« 

^ ^ikxX+ikyy+ika,z~~i J W{C)dC* (42) 



where we denote o;^ = fc^ + A;| + A;| with ki and /c2 given in (p6D, ks = k^ and where W{t) 
stands for an adiabatic series in powers of the time-dependent frequency uj{t) of the modes 
and its derivatives. Up to adiabatic order four (i.e. up to terms involving four derivatives 
of uj{t)) W{t) it is given by, 

^^,/ X Bo Aa Ba Ca Da , . , 

Wit) = cu + ^ + ^ + ^ + ^ + ^ + ^, 43 
tu'^ iu^ uj^ ur ijj^ uj^^ 

where, using the notation V = dV/dt*, 
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^ 32' ^ 128 ' 258 ' 2048 ' 

and ... denote the n adiabatic terms in W{t). Up to adiabatic order zero it is simply 

W(t) = u(t). Observe the two following facts: 

(i) Near the singularity cr = we have W{t) ~ uj(t). This is because the higher adiabatic 
corrections vanish at the singularity. In fact, for the case /? = ±1 in (^61), these higher 
adiabatic corrections naturally vanish. Also, for the case (3 ^ ±1 and under the limit 
prescription introduced in (0), they also vanish. 

1 /9 

(ii) In the flat region prior to the arrival of the waves we have W{t) = uj = {kf + A;| + A;|) 

In that case, since /s = 1, we can use (^) to set t* = t, where without loss of generality 



we choose the value t* = to at t = to. Therefore, the mode solutions (^21) in the fiat region 
reduce to, 



^IV _ \ ^ikxX+ikyy+ikzZ—iuit (45) 



{2Tif2k. 

To see how (|45| ) are related to the fiat modes (p!6D, recall first that the new separation 
constant kz is related to the original k± by the ordinary null momentum relations, i.e., 

a) = A;+ + k^ = k+ — k^. (46) 



In fact, fiat modes (pl| ) have the same formal expression as (^5]) where instead of t there is a 
time-like coordinate T = v+u and instead of z there is a space-like coordinate Z = v—u, and 
coordinates {t, z) and (T, Z) are simply related by a Lorentz boost such that, t = '-/{T+jSZ), 
z = -f{Z + pT), with 

7 = ^. P-^. (47) 

Unfortunately, we use the same notation /? in (^7|) and in (^7|), but note that there is not 
possible confusion. 

It is important to understand, however, that we are constructing a set of mode solutions 
as an adiabatic series in terms of derivatives of the frequency u{t) in the differential equation 



(|39[) . This procedure is similar but not equivalent to the construction of an adiabatic vacuum 
state where the field modes are expanded as an adiabatic series in terms of the derivatives 
of the metric coefficients (see for example |^8| for details). In fact, observe for instance that 



the term fo{t) in ( P^ involves two derivatives of the metric since it is directly related to the 
curvature scalar. Thus, it would be a second order term for an eventual adiabatic vacuum, 
but it is simply a zeroth order term in our adiabatic series in derivatives of uj{t). 



4 Hadamard function in the interaction region 

The key ingredient to calculate the vacuum expectation value of the stress-energy tensor is 
the Hadamard function G^^\x,x'), which is defined as the vacuum expectation value of the 
anticommutator of the field, i.e.. 
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(48) 



This Hadamard function contains non-physical divergence terms which can be subtracted 
by the following point splitting prescription, 



(49) 



where S{x, x') is the midpoint expansion of a locally constructed quantity commonly referred 
as a Hadamard elementary solution (see for example [|^) and given by 



5'(x, x') 



2 
a 



a 



a 



a^i^ ln(/i ^0") 



1 (0) 
-Co cr 

2 ^ . 



ln(yU ^cr) — ^6^2°'''^ 



(50) 



written in Appendix A. We use the 



where the coefficients A'-^-'^p, A*-"^-* 
standard definition for the geodetic hiscalar a{x,x') = {l/2)s'^{x,x'), being s{x,x') the 
proper distance between the points x and x' on a non-null geodesic connecting them. Also, 
(T^(x, x') = {d/dxP-)a{x, x') is a geodesic tangent vector at the point x with modulus s(x, x'), 
being x the midpoint between x and x' on the geodesic. The parameter /i in the logarith- 
mic term of (|50D is an arbitrary length parameter, which is related to the two-parameter 



ambiguity of the point-splitting regularization scheme [^. Then we can compute {T^J) by 
means of the following differential operation. 



(T^,(x)) = lim, D,,,G«(x,x'), 



(51) 



where is a nonlocal differential operator, which in the conformal coupling case (^ = 1/6) 
is given by. 



V 
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(52) 



However, the above differential operation and its limit have no immediate covariant meaning 
because G'^^\x^x') is not an ordinary function but a hiscalar and the differential operator 
V^i, is nonlocal', thus we need to deal with the nonlocal formalism of bitensors (see, for 



example 0, |5T| or the Appendix B of reference [|33| for a review on this subject). 

The regularization procedure (^Sj), however, fails to give a covariantly conserved stress- 
energy tensor essentially because the locally constructed Hadamard function (|^) is not in 
general symmetric on the endpoints x and x' (i.e. it satisfies the field equation at the point 
X but fails to satisfy it at x') (see [^| for details). Thus, to ensure covariant conservation. 



we must introduce an additional prescription: 
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{T,Ax)) = {T^A^))-^^9,.. (53) 

Note that this last term is responsible for the trace anomaly in the conformal coupling 
case, because even though (T^(x)) has null trace when ^ = 1/6, the trace of {T^y{x)) 

is given by {T^) = — a2'^^(x)/(167r^). The regularization prescription just given in ( ^3]) 



satisfies the well known four Wald's axioms [49, 53, 54, 55|, a set of properties that any 



physically reasonable expectation value of the stress-energy tensor of a quantum field should 
satisfy. There is still an ambiguity in this prescription since two independent conserved local 
curvature terms, which are quadratic in the curvature, can be added to this stress-energy 
tensor. In particular, the /i-parameter ambiguity in (^) is a consequence of this (see |^ 



for details). Such a two-parameter ambiguity, however, cannot be resolved within the limits 
of the semiclassical theory, it may be resolved in a complete quantum theory of gravity |^ 



Note, however, that in some sense this ambiguity does not affect the knowledge of the matter 
distribution because a tensor of this kind belongs properly to the left hand side of Einstein 
equations, i.e. to the geometry rather than to the matter distribution. 



5 Procedure to compute {TnJ) 



/iZ// 

After this preliminary introduction on the point-splitting regularization technique, we may 
proceed to calculate the Hadamard function G^^^ (x, x') in the interaction region for the initial 
vacuum state defined by the modes 0^, (^21). The Hadamard function can be written as. 



G^^\x, x') = Mx) <Pl{x') + c.c. (54) 

k 

Note that solutions (pk contain the function h{t*), which cannot be calculated analytically 
but may be approximated up to any adiabatic order as described in (|43|)-(^^. Thus, we have 
the inherent ambiguity of where to cut the adiabatic series. In fact, this is an asymptotic ex- 
pansion, which has a well established ultraviolet limit but it may have convergence problems 
in the low-energy limit. However, observe from (|37D and (|39|) that since: (i) dt/dt* 0, 



(ii) V{t) = uj^it) is bounded and, under the prescription introduced in (0), is also non-zero 
towards the singularity, then the adiabatic series ( ^3]) reduces to ~ u; near the singularity. 
This means that we could cut the adiabatic series ( |45| ) at order zero if we were interested in 
a calculation near the singularity. However, this is only partially true. In fact, it would be 
true if we were only interested in the particle production problem, which essentially would 
involve the evaluation of a Bogoliubov transformation between two different sets of field 
modes. Unfortunately it is not sufficient for the calculation of the vacuum expectation value 
of the stress-energy tensor. This is because G^^^ calculated with h{t*) at order zero does not 
reproduce the short- distance singular structure of a Hadamard elementary solution (pOD in 
the coincidence limit x — > x'. The smallest adiabatic order for the function h{t*) which we 
need to recover the singular structure of G^^^ is order four, basically because our adiabatic 
construction of the mode solutions is similar (but not equivalent) to an adiabatic vacuum 



state (see for details). 

In the mode sum ( ^4|) we use the shortened notation J2k = Io° dk^/k^ I^oo^^x I^^oo^K 
or equivalently J2k = (-^1-^2)"^ c^^i (^^2 I^oo^^s/'^y where the change of variables 
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(p6D and the usual notation (pG]) have been used. Therefore we have, 



1 ^oo ^oo ^oo dkidk2dk3 



G^^Hx x') = / / / ^ X 

2(27r)3LiL2 J^^J~^J-^ Jw{t)W{t') 



e 



-i ^l*,W{C,)dC+ik:,{x-x')+iky{y-y')+ik^{z-z') _^ ^ ^ ^55^ 



We assume that the points x and x' are connected by a non-null geodesic in such a way 
that they are at the same proper distance e from a third midpoint x. We parametrize the 
geodesic by its proper distance r and with abuse of notation we denote the end points by 
X and x', which should not be confused with the third component of [t, z, x, y). Then we 
expand the integrand function in powers of e and we finally integrate term by term to get 
an expression up to e^. The details of such a tedious calculation can be found in ref. 
The result is, 



G'-'\x,x') = A + a6 + C^^aV + D^^^jaVaV + -^|---2A(2Vp^ (56) 



iTl^ I cr 'a 



L 



where L is a logarithmic term defined as L = 27 -|-ln(cr^-R/2) and 7 is the Euler's constant. 
All the involved coefficients A, b, C^p--- depend on which particular space-time, belonging 
to the Szekeres Class of solutions, we are performing the calculations. 

According to (^), the Hadamard function can be regularized using the elementary 
Hadamard solution (jS^) and finally the regularized expression for G^^\x,x') up to order 



IS, 



1 ^ (0)f _[r.(0)A{2) ^J2) 



(57) 



where L is a bounded logarithmic term given by, L = 2'-f + \n{fi'^ fi being the arbitrary 

length parameter introduced in (0), and where the coefficient B is B = b + 3 a'^^ / {32'n''^) . 



From ( p7[ j we can directly read off the regularized mean square field in the "in" vacuum 
state as (0^) = A/2 — a^^^ L / [IGtt'^) . It is important to remark, however, that the term 
D^0^g in ( p7D appears only as a consequence of our approximate procedure of calculating 
the Hadamard function, i.e. using an adiabatic order four expansion for the initial modes 
in powers of the mode frequency uj{t) and its derivatives. Had we used an exact expression 
for the initial modes (or an adiabatic vacuum state ^^), such a term would not appear. 

Now, to calculate the vacuum expectation value of the stress-energy tensor we have to 
apply the differential operator (^) to {^7\). As we have already pointed out, this is not 
straightforward because we work with nonlocal quantities. Note first that the operator ( ^2]) 
acts on bitensors which depend on the end points x and x', but the expression (|5^) for 
depends on the midpoint x. This means that we need to covariantly expand (|57|) in 
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terms of the endpoints x and x' . Also, the presence of quartic terms in ( pTl) gives, after 
differentiation, path dependent terms which must be conveniently averaged. The details of 
such a calculation may be found for instance in p3| , p^ . Then, in the orthonormal basis 
9i = gll"^ dt, 62 = g]J-^ dz, 63 = g]l^ dx, 64 = gljy dy, using the trace anomaly prescription 
(|55D, we may obtain the expectation values {T^^,) in the conformal coupling case and for 
values + e < t < Q oi coordinate t. For values of t < to? {T^^) = 0, and to be consistent 
with the approximation we have used for the space-time geometry, we should require that 
the value of {T^^u) goes smoothly to zero as t ^ tg- In fact, this can be achieved using an 
adequate matching of the line element (|3^) with the flat line element through the interval 
to<t<to + e. 

As a consequence of the logarithmic term in the Hadamard function (|57|) , it will appear 
a similar term in the stress-energy tensor. The argument of this logarithm depends on 
the curvature scalar and thus it will grow unbounded as the flat region is approached. 
However, the coefficient that will appear in front of such a logarithm, depends only on 
locally constructed curvature terms (as can be seen from (|57|)). Therefore, with an adequate 
matching of the space-time geometry, this coefficient will also smoothly vanish towards the 
fiat space region, below t = to- The details of such a matching, however, will not affect the 
main features of the stress-energy tensor, particularly when the singularity is approached. 

Recall that the components of (Tfj^^y) are expressed in the orthonormal thetrad basis 
(6'i, 6'2, 6*3, 6*4) defined upon coordinates (t,z,x,y), with t and z defined in (p9D. These coor- 
dinates allow us to describe symmetrically the head on plane wave collision. Thus, written 
in these coordinates, the stress energy tensor is diagonal. However, an observer at rest at 
the collision center is better described by means of the comoving coordinates T = v + u, 
Z = v — u, such that Z = for the observer's worldline. Observe that the coordinate systems 
(t, z) and (T, Z) are related by the Lorentz transformation (^7|). Therefore, the stress-energy 
tensor in the observer's orthonormal comoving frame will be determined by expressing (T^;y) 

in the orthonormal thetrad basis 61 = g^^ dT, 62 = g^Jz '^'^1 ^3 = 9]Jx '^^1 ^4 = 9yy '^Vi 
defined upon coordinates {T,Z,x,y). This can be directly achieved by transforming {T^J), 
expressed in the basis (6*1, ^^2, ^3, ^'4), by the Lorentz transformation (^T]), i.e.. 



(Tig) = y/5((Tn) + (T22)). (58) 

Observe that in comoving coordinates (T, Z) the collision is not symmetrical anymore and 
the stress-energy tensor acquires a non-diagonal term (TJ2), which indicates the fiux of 
momentum in the Z direction. 



6 Conclusions 

We have given a procedure to calculate the expectation value of the stress-energy tensor 
of a massless scalar field in a family of space-times representing the head on collision of 
two gravitational plane waves throughout the causal past of the collision center and in the 
field state which corresponds to the physical vacuum state before the collision. We have 
considered this particular region essentially because, as we pointed out in a previous work 
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34| , p5| , we could introduce a suitable approximation to the space-time metric (see Fig. 3) 
which not only allow us to dramatically simplify the calculations but also to keep unchanged 
the main physical features, in particular the behaviour of the stress-energy tensor near the 
singularity of the interaction region. 

Unfortunately, a general expression for the stress-energy tensor has not been found. The 
reason is because the approximation (^3]) that we have used in the line element, through- 
out the causal past of the collision center, depends on which particular solution, from the 
Szekeres Class of solutions, we are considering. Nevertheless, one could expect that near 
the singularity the main physical effects would be related to the exponents of the singular 
a{t) terms in the line element (pSl). Thus, the main contribution to the stress-energy tensor 
would be essentially related to these exponents and only trivially sensitive to the remaining 
smooth functions. 
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A Some midpoint expansions 

The coefficients for the midpoint expansion of the locally constructed Hadamard function 
)0f) are: 
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Figure captions 

Fig. 1 The colliding plane wave space-time consists of two approaching waves, regions II and 
III, in a fiat background, region IV, and an interaction region, region I. The two waves move 
in the direction of two null coordinates u and v. The four space-time regions are separated 
by the two null wave fronts u — and v — 0. The boundary between regions I and II is 
{0 < u < 1, V — 0}, the boundary between regions I and III is = 0, < v < 1}, and 
the boundary of regions II and III with region IV is E = < 0, v — 0} U {u = 0, v < 0}. 
Region I meets region IV only at the surface u = v = 0. The singularity in the region 
I corresponds to the hypersurface u"'^ -\- v"^ = 1 and plane wave regions II and III meet 
such a singularity only at V — {u — 1, v — 0} and V' — {u — 0, v — 1} respectively. 
The hypersurfaces u — 1 in region II and ^; = 1 in region III are a type of topological 
singularities commonly referred as folding singularities and they must be identified with V 
and V' respectively. 

Fig. 2 The subset of Cauchy data which affects the evolution of the quantum field along 
the center u — v oi the plane wave collision lies on the segments I^i — {0 < u < s, v — 
0} U {m = 0, <v < s}, where -|- s"^ = 1. Region S is the causal future of this Cauchy 
data (or equivalently, the causal past of the collision center). 

Fig. 3 We change the mode propagation problem for the plane wave collision, in the causal 
past of the collision center, region 5, by a much simpler Schrodinger-type problem which 
consists in: i) performing a Lorentz transformation in order to work with an adequate set 
of coordinates {t,z), ii) eliminating the dependence of the field equation on coordinate z 
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by taking z — and iii) substituting the Cauchy data, which come from the single plane 
wave regions, on segments Ej = {0 < -u < s, v = 0} U {u = 0, < f < s} by much 
simpler Minkowski Cauchy data. This procedure is essentially equivalent to: i) writing the 
line element in coordinates {t,z), ii) modifying the space-time geometry in the causal past 
of the collision center by eliminating the dependence on coordinate z in the line element, 
setting z = and iii) smoothly matching this line element, through plane wave regions II 
and III, with the flat spacetime below the segment {T — 0, —s < Z < s}, where we denote 
T = u + V, Z = V — u. 
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